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1 Preface 



These lecture notes are based on invited lectures at The Second Tah Poe School on Cosmology 
"Modern Cosmology", Naresuan University, Phitsanulok, Thailand, April 17 -25, 2003. See the 
web page 

http : //www . tech . port . ac . uk/staf f web/ gumj udpb/TPC2FIRSTpage . html 
for the details of this school and conference. 

2 Introduction 

The proposal of General Relativity by Einstein in 1915 made it possible to discuss the structure 
of spacetime and the evolution of the universe in terms of physical laws. In 1922 Friedmann 
found the existence of expanding/collapsing cosmological solutions by solving the Einstein field 
equations. In 1929 Hubble discovered the expansion of the universe by the observations in the 
redshift of galaxies, as the Einstein theory predicts. In 1946 Gamov and his collaborators showed 
that the universe must begin in a very hot and dense state from the theory of nucleosynthesis. 
They also predicted that the present universe should be filled with the microwaves with black 
body radiations. In 1965 Penzias and Wilson discovered microwave background radiations which 
coincide well with theoretical predictions by Gamov et ai. These strong observational evidences 
have made people believe that the universe started out from a hot and dense state, which is 
called as the standard big-bang model. 

In standard big-bang cosmology the state of the universe is characterized by the radiation- 
dominated or the matter-dominated stage. These correspond to the decelerated expansion of the 
universe where the second derivative of the scale factor is negative. Meanwhile this decelerated 
expansion is not sufficient to solve a number of cosmological problems such as flatness and 
horizon problems which plagues in the standard big-bang scenario^. In order to overcome these 
fundamental problems, it is required to consider an epoch of accelerated expansion in the early 
universe, i.e., inflation. 

The basic ideas of inflation were originally proposed by Guth P and Sato ^ independently 
in 1981, which is now termed as old inflation^. This corresponds to the de-Sitter inflation which 
makes use of the first-order transition to true vacuum. However, it has a serious shortcoming 
that the universe becomes inhomogeneous by the bubble collision soon after the inflation ends. 
The revised version was proposed by Linde and Albrecht and Steinhardt K in 1982, which 
is dubbed as new inflation. This corresponds to the slow-roll inflation with the second-order 
transition to true vacuum. Unfortunately this scenario also suffers from a flne-tuning problem of 
spending enough time in false vacuum to lead to sufficient amount of infiation. In 1983 Linde 
considered the variant version of the slow-roll inflation called chaotic inflation, in which initial 
conditions of scalar fields are chaotic. According to this model, our homogeneous and isotropic 
universe may be produced in the regions where inflation occurs sufficiently. While old and 
new infiation models are based on the assumption that the universe was in a state of thermal 
equilibrium from the beginning, chaotic inflation can occur even without such an assumption. 
In addition chaotic inflation can start out in the regime close to Planck density, thereby solving 
the problem of initial conditions. 

^In the next section I will explain how inflation solves these problems. 

^Note that a specific version of the inflationary scenario-so called inflation- was proposed by Starobinsky 
one year earlier Nevertheless this did not explicitly point out the virtue of inflation as in the paper of Guth. 
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Many kinds of inflationary models have been constructed in these twenty years (see e.g., 
ref. [3 IB])- In particular, the recent trend is to construct consistent models of inflation based on 
superstring or supergravity models. See ref. for a review in such a direction. 

The inflationary paradigm not only provides an excellent way in solving flatness and hori- 
zon problems but also generates density perturbations as seeds for large scale structure in the 
universe In fact inflation provides a causal mechanism to generate nearly scale- invariant 
spectra of cosmological perturbations. Quantum fluctuations of the field driving inflation-called 
inflaton- are typically frozen by the accelerating expansion when the scales of fluctuations leave 
the Hubble radius. Long after the inflation ends, the scales cross inside the Hubble radius 
again. Thus perturbations imprinted during inflation can be the origin of large-scale structure 
in the universe. In fact temperature anisotropics observed by the COBE satellite in 1992 exhibit 
nearly scale-invariant spectra as predicted by the inflationary paradigm. Recent observations of 
WMAP also show strong evidence for inflation 



3 The standard big-bang cosmology and its problems 

3.1 The standard big-bang cosmology 

The standard big-bang cosmology is based on the cosmological principle [H], which means that the 
universe is homogeneous and isotropic on large distance. Then the metric takes the Friedmann- 
Robertson- Walker (FRW) form: 

Here a[t) is the scale factor with t being the cosmic time. The constant K is the spatial cur- 
vature, where positive, zero, and negative values correspond to closed, flat, and open universes, 
respectively. 

The evolution of the universe is dependent on the material within it. This is characterized 
by the equation of state between the energy density p{t) and the pressure p{t). Typical examples 
are : 

p = /c/S, radiation, (2) 
p = 0, dust. (3) 

In order to know the dynamical evolution of the universe, it is required to solve the Einstein 
equations in General Relativity. The Einstein equations are expressed as ^1 IT^ 

G^,y = R^u -'^g^,yR = '^TiGT^v - Kg^v, (4) 

where i?, T^y, and G are the Ricci tensor, Ricci scalar, energy momentum tensor, gravi- 
tational constant, respectively. The Planck energy, m^i = 1.2211 x 10^^ GeV, is related with G 
through the relation nipx = /G)^^"^ . Here fi and c are the Planck's constant and the speed of 
light, respectively. Hereafter we use the units h = c = 1. A is a cosmological constant originally 
introduced by Einstein. 



(1) 
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For the background metric with a neghgible cosmological constant, the Einstein equations 
yield 



p + 3H{p + p) = 0, (6) 

where a dot denotes the derivative with respect to t, and H = a/a is the Hubble expansion rate. 
Eqs. and © are so called the Friedmann and Fluid equations, respectively. Combining these 
relations gives the following equation, 

a 4:TT , 

The Friedmann equation @ can be rewritten as 

fi-l = ;^. (8) 

where 

n = ^, with p, = ^-^. (9) 
Pc on 

Here the density parameter is the ratio of the energy density to the critical density. When 
the spatial geometry is flat {K = 0), the solutions for Eqs. Q and © are 

Radiation dominant : 

Dust dominant : a oc p (X a . (11) 

In these simple cases, the universe is expanding deceleratedly (a < 0) as confirmed by Eq. 

3.2 Problems of the standard big-bang cosmology 
3.2.1 Flatness problem 

In the standard big-bang theory with a < 0, the a^if^(= d^) term in Eq. Q always decreases. 
This indicates that tends to shift away from unity with the expansion of the universe. However, 
since present observations suggest that Q is within the order of magnitude of one , O needs 
to be very close to one in the past. For example, we require — 1| < O{10~^^) at the epoch of 
nucleosynthesis |H1II2 and < ©(lO"*^^) at the Planck epoch This is an extreme fine- 

tuning of initial conditions. Unless initial conditions are chosen very accurately, the universe 
soon collapses, or expands quickly before the structure can be formed. This is so called the 
flatness problem. 
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3.2.2 Horizon problem 



Consider a comoving wavelength, A, and also a physical wavelength, oA, which is inside the 
Hubble radius, {i.e., aX < H^^). The standard big-bang cosmology is characterized by the 
cosmic evolution of a cx with < p < 1. In this case the physical wavelength grows as aA oc t^, 
whereas the Hubble radius evolves as H^^ oc t. Therefore the physical wavelength becomes much 
smaller than the Hubble radius with the passage of time. This means that the region where the 
causality works eventually becomes the only small fraction of the Hubble radius. 

To be more precise, let us first define the particle horizon Dnit) where the light travels from 
the beginning of the universe, t = t^, 



Here dnit) corresponds to the comoving distance. Setting = 0, we find dnit) = 3t in 
the matter-dominant era. We observe photons in the Cosmic Microwave Background (CMB) 
which are emitted at the time of decoupling. The particle horizon at decoupling, Dnitdec) = 
o(^dec)'^H(^dec)> Corresponds to the region where photons could have contacted causally at that 
time. The ratio of dnitdec) to the particle horizon today, dnito), where to is the present time, 
is estimated as 



This result implies that the causality regions of photons are restricted to be small. In fact the 
surface of the last scattering surface only corresponds to the angle of order 1°. Observationally, 
however, we see photons which thermalize to the same temperature in all regions in the CMB 
sky. This is so called the horizon problem. 

3.2.3 The origin of large-scale structure in the universe 

The COBE satellite observes the anisotropics in the last scattering surface, whose amplitudes 
are small and close to scale-invariant. These fiuctuations spread so large a scale that it is 
practically impossible to generate them between the big bang and the time of the last scattering 
in the standard cosmology. This problem is almost equivalent to the horizon problem mentioned 
above, i.e., the standard cosmology can not provide satisfactory explanation for the origin of 
large-scale structure. 

3.2.4 Monopole problem 

According to the view point of particle physics, the breaking of supersymmetry leads to the 
production of many unwanted relics such as monopoles, cosmic strings, and topological defects 
[7]. The string theories also predict super symmetric particles such as gravitinos, Kaluza-Klein 
particles, and moduli fields. 

If these particles exist in the early stage of the universe, the energy densities of them decrease 
as a matter component (~ a~^). Since the radiation energy density decreases as ~ in the 
radiation-dominant era, these massive relics could be the dominant materials in the universe, 
which contradicts with observations. This problem is generally called as the monopole problem. 




with 




(12) 




(13) 
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4 Idea of inflationary cosmology 



The problems in the big-bang cosmology lie in the fact that the universe always exhibits the 
decelerating expansion. Let us assume the existence of a stage in the early universe with an 
accelerated expansion of the universe, i.e., 

a>0. (14) 

From the relation ((Jj) this corresponds to the condition 

p + 3p<0. (15) 

The condition (|14() essentially means that d (= aH) increases during inflation. Then the comov- 
ing Hubble radius, {aH)~^, decreases in the inflationary phase. This property is the key point 
to solve the cosmological puzzles in the standard big-bang cosmology as I will show below. 



4.1 Flatness problem 

Since the a?H^ term in Eq. increases during inflation, Q rapidly approaches unity. After the 
inflationary period ends, the evolution of the universe is followed by the conventional big-bang 
phase and — 1| begins to increase. In spite of this, as long as the inflationary expansion occurs 
sufficiently and makes Q, very close to one, Vt stays of order unity even in the present epoch. 



4.2 Horizon problem 

Since the scale factor evolves as a oc with p > 1 during inflation, the physical wavelength, aA, 
grows faster than the Hubble radius, H^^{(x t). Therefore the physical wavelength is pushed 
outside the Hubble radius during inflation. This means that the region where the causality works 
is stretched on scales much larger than the Hubble radius, thus solving the horizon problem. 

Of course the Hubble radius begins to grow faster than the physical wavelength after inflation 
(radiation/matter dominant era). In order to solve the horizon problem, it is required that the 
following condition is satisfied for the comoving particle horizon: 

» / -TIT • (16) 



This implies that the comoving distance that photons can travel before decoupling needs to be 
much larger than that after the decoupling. According to the precise calculation, it is achieved 
when the universe expands about e^^ times during inflation jSl I15j. 



4.3 The origin of the large scale structure 

The fact that the comoving Hubble radius decreases during inflation makes it possible to generate 
the nearly scale-invariant density perturbations on large scales. Since the scales of perturbations 
are within the Hubble radius in the early stage of inflation, causal physics works to generate 
small quantum fluctuations. After a scale is pushed outside the Hubble radius (i.e., the first 
horizon crossing) during inflation, the perturbations can be described by classical ones. When 
the inflationary period ends, the evolution of the universe is followed by the standard big-bang 
cosmology, and the comoving Hubble radius begins to increase. Then the scales of perturbations 
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cross inside the Hubble radius again (the second horizon crossing), after which causahty works. 
The small perturbations imprinted during inflation appear as large-scale perturbations after 
the second horizon crossing. In this way the inflationary paradigm naturally provides a causal 
mechanism in generating the seeds of density perturbations observed in the CMB anisotropics. 

4.4 Monopole problem 

During the inflationary phase (p + < 0), the energy density of the universe decreases very 
slowly. For example, when the universe evolves as a oc with p > 1, we have H (x oc a~^^^ 
and p oc ot'^I'^ . Meanwhile the energy density of massive particles decreases much faster (~ a~'^), 
these particles are red-shifted away during inflation, thereby solving the monopole problem. 

Of course, we have to worry for the case where these unwanted particles are produced after 
inflation. In the process of reheating followed by inflation, the energy of the universe can be 
transferred to radiation or other light particles. At this stage unwanted particles must not be 
overproduced in order not to violate the success of the standard cosmology such as nucleosyn- 
thesis. Generally if the reheating temperature is sufficiently low, the thermal production of 
unwanted relics such as gravitinos can be avoided. 

5 Inflationary dynamics 

The scalar fields are the important ingredients in particle physics theories. Consider a homoge- 
neous scalar field, (/>, called inflaton, whose potential energy leads to the exponential expansion 
of the universe. The energy density and the pressure density of the infiaton can be described, 
respectively, as 

p = \'P^ + Vi<P), p=^^^-Vi^), (17) 

where V{(p) is the potential of the infiaton. Substituting Eq. ifTTjl for Eqs. © and ©, we get 

Svr 



(18) 



4> + 3H<j) + V'{^) = 0, (19) 

where = SvrG = Svr/mpj, and we neglected the curvature term K"^ /a? in Eq. H18|) . 

During inflation, the relation ()15() yields (jP' < V{(p), which indicates that the potential energy 
of the inflaton dominates over the kinetic energy of it. Therefore a flat potential of the inflaton 
is required in order to lead to sufficient amount of infiation. Imposing the slow-roll conditions: 
^<p^ <C V{(1)) and (f> ^ 3H<p, Eqs. (|18|) and (fT!l|) are approximately given as 



^ i^V{^), (20) 
3m 



pi 



3H<j) ~ -y'(0). (21) 
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Defining the so-called slow-roll parameters 




rj = 



m-l^ V" 



(22) 



we can easily verify that the above slow-roll approximations are valid when 



e<Cl, |r/|<l. 



(23) 



The inflationary phase ends when e and \r]\ grow of order unity. A useful quantity to describe 
the amount of inflation is the number of e-foldings, defined by 



where the subscripts i and / denote the quantities at the beginning and the end of the inflation, 
respectively. 

In order to solve the flatness problem, Q is required to be — 1| < 10^^*^ right after the 
end of inflation. Meanwhile the ratio |0 — 1| between the initial and final phase of inflation is 
given by 



where we used the fact that H is nearly constant during inflation. Assuming that |0j — 1| is of 
order unity, the number of e-foldings is required to be A'^ > 70 to solve the flatness problem. We 
have similar number of e-foldings from the requirement to solve the horizon problem. 

6 Models of inflation 

So far we did not mention the form of the inflaton potential, V{(p). Originating from the old 
inflationary scenario proposed by Guth 1 and Sato [2], we now have varieties of inflationary 
models : R^, new, chaotic, extended, power-law, hybrid, natural, supernatural, extranatural, 
eternal, D-term, F-term, brane, oscillating, trace-anomaly driven,..., etc. 

These different kinds of models can be roughly classified in the following way 16 . The first 
class (type I) is the "large field" model, in which the initial value of the inflaton is large and 
it rolls down toward the potential minimum at smaller cf). Chaotic inflation |Sj is one of the 
representative models of this class. The second class (type II) is the "small field" model, in 
which the inflaton field is small initially and slowly evolves toward the potential minimum at 
larger (f). New inflation and natural inflation |17j are the examples of this type. In the 

first model, the second derivative of the potential V^'^^ {(f)) usually takes positive values, whereas 
can change its sign in the second model. The third one (type III) is the hybrid (double) 
inflation model ^Sj, in which inflation typically ends by the phase transition triggered by the 
presence of the second scalar field (or by the second phase of inflation after the phase transition) . 

Let us briefly recap each type of inflationary models. 




(24) 




(25) 
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6.0.1 Chaotic inflation— an example of type I 

The chaotic inflation model is described by the quadratic or quartic inflaton potential, 

y(<A) = ^mV, or V{cl>) = -^M^. (26) 

The "chaotic" means that initial conditions of inflaton are distributed chaotically. According 
to this scenario, the region which undergoes the sufficient amount of inflation gives rise to our 
universe. 

Let us investigate the evolution of the universe in the case of the quadratic potential. Then 
Eqs. (EOI) and ^ read 

^2^47rmV^ 2,H^ + m^^Q. (27) 

Combining these relations gives the following solutions: 

"^"^-{t, (28) 



2V37r 



a ~ Oj exp 



3 mpi V 4V37r 



(29) 



where (j)i is an integration constant corresponding to the initial values of the inflaton. We 
find from the relation H29|) that the universe expands exponentially during the initial stage of 
inflation. With the increase of the second term in the square bracket of Eq. (|29j) . the expansion 
rate slows down. Since the slow-roll parameters are expressed as 



m 



V 



pi 



(30) 



the inflationary period ends around |(/>| ~ mpi/v47r, after which the system enters a reheating 
stage. The total amount of inflation is approximately expressed as 



N 



- -) 



(31) 



In order to lead to sufficient inflation > 70, we require the initial value to be (f)i > SrUpi. 
The inflaton mass, m, is constrained by the amplitude of density perturbations observed by the 
COBE satellite. In order to fit observations, m is required to be jH] 



m 



10"^mpi 



(32) 



In the quartic potential case, the self-coupling is constrained to be A ~ 10 by the similar 
arg ument [ZIIH]- 
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Figure 1: The schematic illustration of the potential of the chaotic inflation model. This belongs 
to the class of the "large field" model. 



6.0.2 Natural inflation an example of type II 

Natural inflation model is characterized by Pseudo Nambu-Goldstone bosons (PNGBs) 
which appear when an approximate global symmetry is spontaneously broken. The PNGB 
potential is expressed as 



m 



1 + cos I y 



(33) 



where two mass scales m and / characterize the height and width of the potential, respectively. 
The typical mass scales are of order / ~ rripi ~ 10^^ GeV and m ~ mQUT ~ 10^^ GeV, which 
are expected by particle physics. 

Let us consider the case where the inflaton is initially located in the region, < (j) < vr/, 
and inflation occurs while the inflaton slowly evolves toward the potential minimum at (j) = vr/. 
The slow-roll parameters are 



m 



pi 



167r/2 



sin(0//) 



1 + cos((/>//) 



m^i cos{(t>/f) 
'87r/2l + cos(0//)' 



(34) 



Note that e and ry depend on /, but not on m. Inflation starts out from the regime where (f) is 
close to zero. The system enters a reheating stage when the inflaton begins to oscillate around 
4> = tt/. 

One typical property in the type II model is that the second derivative of the inflaton 
potential can change its sign. In the natural inflation, V^'^\(l)) is negative when inflaton evolves 
in the region < (p < vr//2. This leads to the enhancement of inflaton fluctuations by 
spinodal (tachyonic) instability [W\ \2()\ I21j . When the particle creation by spinodal instability 
is neglected, the number of e- foldings is expressed by PT| 



N 



In 



m: 



pl 



sm 



sin(0,/2/) 



(35) 
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where (pi and (pf are the initial and final values of the inflaton during inflation, respectively. In 
order to obtain the sufficient amount of e-foldings satisfying N > 70, the initial value of the 
inflaton is required to be (j){ti) < O.lmpi for the mass scale / ~ mpi. 

6.0.3 Hybrid (Double) inflation an example of type III 

The inflationary model building in the presence of multiple scalar fields is the recent trend from 
the viewpoint of particle physics 0. Let us consider the Linde's hybrid inflation model |18j . 
described by 

V = 2(x'-^J + + \rn'cP' ■ (36) 

When (/)^ is large, the field tends to roll down toward the potential minimum at x = 0. In this 
case, the potential is effectively described by a single field, 

- ^ + lrn'<P' ■ (37) 

Inflation does not end for the potential ()37() . In the presence of the x field, however, the mass 
of the field x becomes negative for (p < (p^ = M/g. Then the field begins to roll down to 
one of the true minima at = and x = ±M/\/A (see Fig. The original version of the 
hybrid inflation ^HI corresponds to the case where the inflation soon comes to an end after the 
symmetry breaking {(p < (pc) due to the rapid rolling of the fleld x- this case the number of 
e-foldings acquired in double inflation can be approximately estimated by using the potential 
(EZI): 

2ttM^ 



iV~-^^ln^, (38) 
Am^rripj (pc 
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/ 



Figure 3: The schematic ihustration of the potential of the hybrid (double) inflation model. 
This model is characterized by multiple scalar fields. 

where is the initial value of inflaton. 

If the condition, » AM^, is satisfied, the mass of the field x is "light" relative to the 
Hubble rate around (p = (j)c, thereby leading to the second stage of inflation for (p < (j)^ . This 
corresponds to the double inflationary scenario. 

7 Density perturbations from inflation 

In this section we study cosmological perturbations generated from inflation. For simplicity we 
shall consider the single-field model of inflation. The most general form of the line element that 
describes scalar metric perturbations is written as [231 124j 



where t = J a^^dt is the conformal time. The scalar quantities A, ijj, B, and E are the functions 
of space and time coordinates. We denote the three-dimensional covariant derivative of a function 
/ with respect to some coordinate i by /j. In the linear perturbation theory of cosmological 
perturbations, there exist some redundant degress of freedom in the equations of motions which 
can be fixed without affecting the physics. It is important to choose a proper gauge condition 
so that it allows for an simplified mathematical analysis and an easier physical interpretation 
(see e.g., Ref. [^S] for details). Among several gauge conditions, the longitudinal gauge defined 
by the gauge condition B = E = is often used : 



ds'^ = a2(r){-(l + 2A)dT'^ + 2B,i dx'dr + [(1 - 2i;)jij + 2E,ij ]dx'dx^} , 



(39) 



ds 



,2 



a2(r)[-(l + 2^)dT^ + (1 - 2^)jijdx'dx^] 
-(1 + 2^)dt^ + a2(t)(l - 2'^)5ijdx'dx^ . 



(40) 
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The potentials ^ and ^' can be constructed as gauge-invariant variables, whose forms are invari- 
ant under infinitesimal coordinate transformations (see Ref. j24j for details). When the spatial 
part of the energy- momentum tensor is diagonal, 5Tj ~ (^j, ^ coincides with ^. 
Let us consider the linearized Einstein equation with A = : 



6G^^ = SttGSTi;. 

The perturbed Einstein tensor can be evaluated in the longitudinal gauge as follows 



(41) 



5G° 



6G] 



2a 



-3W(W$ + ^') + 



-2a' 



1. 



(42) 



where 7i = a' /a and D = ^ — ^, with a prime being a derivative with respect to r. In the 
case where the universe is filled with a minimally coupled scalar field, </>, the energy momentum 
tensor takes the following form 



T,^ = -(V0)2 + 



^(V<^)2 + y(^) 



We decompose the field (j) iiito homogenous and fluctuational parts as (j) = (j)Q + c 
energy momentum tensor is decomposed into background and perturbed parts as 



where 



and 



=(0) j^M + ^T^M ^ 



-^<A? + y(<Ao) 



6]. 



(43) 
Then the 

(44) 
(45) 



5T^ = 




</.' + V'{(l))a'^5cb 












st; = 









(46) 



Note that the relation JTj ~ (5j holds in this case. This leads to the relation Z) = ^> — ^' = Oby 
considering the ij part in Eq. H41|). Then the perturbed Einstein equation H41|) can be written 

as 



- m^' - {W + 2^2)$ = y (0o(5(/. + V'{(i))a^ 



5(1) - V\(l))a^5(l)) , 



(47) 
(48) 
(49) 
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where k? = SttG. We used the background equations in deriving these equations. Combining 
Eqs. (|47j) - (|49)) gives the equation for the gauge-invariant scalar field: 

6(1}" + 2n6(p' - V'^6(p + y"((/>)a2(50 = A^'q^' - 2V'{(l))a'^^. (50) 

Note that the effect of gravitational perturbations appears in the rhs of Eq. (|5U|) . 

Eliminating the scalar field perturbation in Eqs. H47j) - H49j) . we find the following equation for 
gauge-invariant cosmological perturbations: 

+ 2(a/(/>'o)'(a/0[))-i$' - + 2cl)'o{n/ <!)')' ^ = . (51) 

Introducing new variables 

u = (a/(/>'o)$ , z = n/acj)'o , (52) 

Eq. reduces to 

u" - V^u - {z"/z)u = . (53) 
Making Fourier transformations, the mode function Uk with a wave number k satisfies^ 

u'l + (^2 _ z"/z)uk = . (54) 
In the long-wavelength limit, <^ z" /z, we find the analytic solution 

u = z(^ci + C2 J , (55) 

where ci and C2 are integration constants. Making use of Eqs. (|48|) . (|52j) . and (|55|) . we have 

<^ = -^(ci- 2c2 J adtj + 2c2 , (56) 



— ci - 2c2 / adt . (57) 



a 

The solution H56|) indicates that there exists a gauge-invariant conserved quantity (comoving 
curvature perturbation) in the long-wavelength limit: 

'^^^'Tfr^Hj ^^^^ 

which was firstly introduced by Bardeen |26j . 

During inflation, the perturbations with a comoving wavenumber k are pushed outside the 
Hubble radius by the accelerated expansion of the universe. This makes the wavenumber toward 
the long wavelength regions: k'^ <C z" /z, in which case we have the analytic solutions of Eqs. (|56|) 

^Note that the comoving wavelength A is related with k by the relation \ — 2Tv/k. 



14 



and (|57() . Since the terms which include the ci coefficient rapidly decrease and the relation 
^ holds during the inflationary phase, Eqs. (|56jl and l|57jl yield 



2c2 ( — / adt 



2C2 



H' 



a{H-^ydt 



(59) 



Sep — 2c2(j)oa ^ I adt ~ 2c2(t)o/H . 
Combining Eqs. and ((HOJ, we have 



7^ 



(60) 



(61) 



where we introduced the subscript "A;" for the fluctuation, 5(j). The H5(j)k/4'o term should be 
evaluated when the wavelength of the perturbation crosses the Hubble radius: /c^ = z" /z. Once 
the scales are pushed outside the Hubble radius during inflation, the evolution of perturbations 
can be described by Eqs. and ((57j) . Long after the inflation ends, the scale of the perturba- 
tions reenters the Hubble radius again, which generates the seeds for the galaxies and clusters. 
Thus the inflationary paradigm naturally explains the origin of large-scale structure. 
The spectrum of the comoving curvature perturbation is defined as 



2vr2 ^2 



(62) 



From Eqs. (|59|) and (|6Ujl we have $ ~ eH5(j)k/4>o with e = —H/H^ being the slow-roll parameter 
defined in eq. ((211). Taking note that |4</)'o^>'| <C \2V'{(l))a^<^\ on super-Hubble scales and making 
use of the slow-roll relation 3H(]) ~ —V'{<j)), Eq. (|5U|) can be approximately written as 



+ 2n6(P'k + + aV 



Qa^eH^ 



0. 



Redefining a field perturbation, 6xk = a ^6(j)k, this equation reduces to 



^Xk + 



^Xk = with z^^ = — + 9e — 3// . 



(63) 



(64) 



where we used the fact that a" /a = a^{2- e)H^ ~ (2 + 3e)/T^. When is a real number, the 
solution for this equation can be written in terms of the Hankel functions of the first and second 
kind: 



^Xk 



T Cl 



{k)Hi'\-kT)+C2{k)Hi^H-kr) . 



(65) 



Quantum fluctuations are generated in the ultraviolet regime {k ^ aH), which is required 
to take the plane-wave form e~'^^'^/V2k for —kr ^ 1. Matching Eq. H65|) with this quantum 
fluctuation, we have ci{k) = ^^e^^^^^)^ and C2(A;) = 0, thereby yielding 



^Xk 



THl^\-kT). 



(66) 
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Since the Hankel function takes the form Hy{-kT < 1) ~ ^;^e"*2 2''-2 (r(i/)/r(3/2))(-A;r) 
on super-Hubble scales, the amphtude of the fluctuation 5(t)k is estimated as 



H / k 
7W \aH 



(67) 



Substituting this relation for eq. H62() . we finally get the spectrum of the curvature perturbations 
on super-Hubble scales: 

Here the spectral index n-jz. is given as 

_ d\nPn „ ^„ 

^^-' = ^h^ = -'^+'^- 

When n-ji = 1, the spectrum of the curvature perturbation is scale-invariant. Since the slow-roll 
parameters e and rj are much smaller than unity [see eq. (|2,''{|) ]. scalar perturbations generated 
in standard inflation are close to scale-invariant. This property is consistent with the recent 
observation of WMAP When n-ji < 1 and nn > 1, we call the red-tilted spectrum and 

the blue-tilted spectrum, respectively. For example the chaotic inflation model exhibits the 
red-tilted spectrum, while the hybrid inflation model shows the blue-tilted spectrum :.8:- This 
property is useful to constrain on inflationary models from observations. 

The spectrum (|68|) corresponds to the adiabatic perturbation generated in the single-field 
inflationary scenario. In the case of the multi-field inflation, this is modified due to the presence 
of isocurvature (entropy) perturbations |22j. In particular the curvature perturbation TZ is 
no longer a conserved quantity in the multi-field context |28j. Therefore we have to estimate 
the curvature perturbation at the end of inflation instead of using its value at the horizon 
crossing. See ref. ^T. for the detailed numerical analysis of correlated adiabatic and isocurvature 
perturbations in concrete models of double inflation. 



8 Reheating after inflation 

After the end of inflation, the universe enters a reheating stage, during which the potential 
energy of the inflaton is transferred to radiation and the universe is thermalized. In the original 
version of the reheating scenario which is now called old reheating, the decay of the inflaton 
is described by the perturbation theory This process is not efficient for the success of 

the GUT scale baryogenesis scenario. In contrast, it was later found that the existence of the 
nonperturbative stage called preheating can lead to the explosive particle production in the early 
stage of reheating [301 1^ ■ 



8.1 Old reheating 

The reheating stage turns on when the inflaton reaches the potential minimum and begins 
to oscillate. Let us first study the background evolution during reheating for the polynomial 
potential 

y(^) = J-m^m^i ( ^] , (70) 
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where n is a positive integer. Making use of the time-averaged relation 



{4>'')T = 2n{V{cf>))T, (71) 
the time evolution of the scale factor and the Hubble parameter are 

a oc , Hoc^-, (72) 

3n t 

where we used eqs. (|18() and ()19p . When n = 1 and n = 2, the universe evolves as matter- 
dominant (a oc t^/^) and radiation-dominant (a oc t^^"^), respectively. 

Let us consider the case of n = 1 in order to understand the basic picture of reheating. In 
this case, the evolution of the inflaton is described by the sinusoidal oscillation with a decreasing 
amplitude ^(t): 

^ = ^{t)smmt, i>(t) = ^£^. (73) 

Then the energy density of the inflaton decreases as 

p = i,/,2 + y(<^) ~ im2$2 Of . (74) 

The simple estimations of Eqs. ()73p and H74|) should be modified in the presence of the infiaton 
decay. Without this decay, the universe will evolve toward emptier and colder states because of 
the redshift of the energy density. In the old reheating scenario ^H] where the single infiaton 
decay is considered, the infiaton <j) is coupled to a scalar field x & fermion field ip with the 
interacting Lagrangian, 

Ant = -o-'/'X^ - hcpijj^ , (75) 

where a and h are coupling constants. Note that we assume that the bare masses of x &iid ip 
are much smaller than the infiaton mass m and neglect them. We include quantum correction 
effects to the equation of infiaton as 



+ 3H(^ + 



+ n(m) = 0, (76) 



where n(m) is the polarization operator of infiaton. Although the real part of n(m) is small 
relative to m?, 11 (m) acquires an imaginary part, 

Im n(m) = —mT , (77) 

where F = T{(j) xx) + ^{4' ~^ V'V') is the total decay rate of infiaton. The perturbative decay 
rates of (p ^ x and 4> ^ are expressed as [2^1 

^<*^«) = 8^- = (78) 

Note that these relations are valid in the case of T ^ m, which implies the relations: cr^ <C 
and /i^ <C 1. We cannot use Eq. ((75|l in the nonperturbative stage of reheating (preheating). 
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Combining Eqs. ()76() and ()77() and setting = <I>(t) exp{imt), we find 



$(t) = -|H^exp(-rt/2). (79) 



The same solution can be obtained by adding the phenomenological decay term T(j) to the 
equation of inflaton: 

4> + 3H,j) + T^ + rr?(\) = . (80) 

The relation, F < 3//, holds in the initial stage as long as coupling constants a and h are small. 
Since the Hubble parameter decreases as H oc the fraction of produced particles to the 
total energy density becomes important when drops less than F. The energy density of the 
universe at this time can be estimated by setting = {3H)'^ = 247r/>/mpj, as 

If we assume that this energy density is rapidly transferred to light particles and newly-produced 
particles acquire the radiation energy with temperature Tr by instant thermalization, we obtain 
the relation 

where g*{ ^ 100) is the effective number of degrees of freedom at T = Tr. Then the reheating 
temperature is estimated as 



rfi<0.1^rmpi. (83) 
Taking into account the relation, T <^ m < lO^^mpi, we find 

Tr « 10~^mpi . (84) 

This value is much smaller than the GUT scale temperature Tgut = 10~'^mpi, which indicates 
that the GUT scale baryogenesis does not work well in the old reheating scenario. 



8.2 Preheating 

In the early 1990s, it was realized that the inflaton decay may have started in a much more 
explosive process, called preheating jHOl 1^ before the perturbative decay (see also refs. [HSl 
ESI EH)- At this stage, the fluctuations of scalar particles can grow quasi-exponentially by 
parametric resonance. The process of reheating consists of three stages: 

(i) Preheating stage where particles are nonperturbatively produced 

(ii) Perturbative stage where the inflaton decay is described by the Born process 

(iii) Thermalization of produced particles 

Among variant models of inflation, chaotic inflation has an advantage that the dynamics of 
preheating is simply analyzed. Kofman, Linde, and Starobinsky investigated the structure of 
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resonance in preheating both in massive and self- interacting inflationary models Hereafter 
we shall review the dynamics of preheating for the two-field model with a massive inflaton, 

y(</.,x) = ^mV + ^5W. (85) 

We assume that the spacetime and the inflaton (p give a classical background and the scalar field 
X is treated as a quantum field on that background. Expanding the scalar field x as 

X = I {akXk{t)e-^''-'' + aixme^"^--) d'k , (86) 

and adopting the flat Friedmann-Robertson- Walker (FRW) metric, the each Fourier component 
Xk{t) satisfies the following equation of motion, 

Xk + SHxk+{^+g^<P^]xk = 0. (87) 



Introducing a new scalar field = a^/^Xfci Eq. (|87j) yields 



Xk + ^k^k = , 



where 



2 fc2 2 ,2 3 /2a 

We neglect the last term in Eq. H89|) for analytic investigations since this term becomes gradually 
negligible during the reheating phase. 

Then Eq. H89|) is reduced to the well known Mathieu equation 

(fx 

-^ + {Ak-2qcos2z)Xk = 0, (90) 

where z = mt and 

Ak = 2q + , (91) 



4m ^ 

The strength of the resonance in Eq. (|9Uj) depends on the variables of Ak and q, which is shown 
in a stability-instability chart of Fig. [l] In the unstable region (the lined region in Fig. Xk 
grows exponentially as Xk oc exp(/Xfc2:) with the Floquet index fik and particles with momentum 
k are produced. For small q (< 1), the width of the instability band is small and few /c-modes 
grow by this resonance. This is called the narrow resonance. On the other hand, for the large 
q{'^ 1), the resonance can occur for a broad range of the momentum A;-space. Since the growth 
rate of the x-particle becomes larger with the increase of the variable q, this resonance provides 
more efficient particle production than in the narrow one. This is called the broad resonance |,Slj . 
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Note that the initial amplitude of the inflaton and the coupling constant g play important roles 
to determine whether the resonance is narrow or broad. Since the inflaton mass is constrained 
to be small as m ~ lO^^rripi by COBE normalizations, the large resonance parameter, q ^ 1, 
can be easily achieved for the coupling g > 10~^ with the initial amplitude, ^{t-i) ~ 0.2mpi. The 
allowed region on the Mathieu chart is determined by Eq. (|91|) as 

Ak>2q. (93) 

Hence the broadest resonance is given by the limit line Aj. = 2q. 

As is found in Eq. (|91j) . particles with low momenta are mainly produced. Considering 
the nonadiabatic condition, duj^/dt ^ ojf,, for the particle production, the maximum comoving 
momentum is roughly estimated as [HHl 




This indicates that large values of g and $ lead to the production of particles with high momenta. 

When q is sufficiently large initially ((? ^> 1), the resonance of each mode occurs stochastically 
[33] . In this case, the frequency uJk decreases by cosmic expansion and ujk drastically changes 
within each oscillation of the inflaton field, so the phases of x held at successive moment of 
= are not correlated each other. In the first stage of preheating, the x fields cross large 
number of instability bands. The periods when they are in instability bands are so short that 
the resonance can not occur efficiently relative to the Minkowski spacetime case. Nevertheless, 
the X fluctuations can still grow quasi-exponentially. As q becomes smaller, cosmic expansion 
slows down, and the fields stay in each resonance band for a longer time. When the variables 
decrease below the lower boundary of the first resonance band by the expansion of the universe, 
particle productions come to an end. 

We have to note that there is another mechanism which terminates the resonance. Taking 
into account the backreaction effect of created x particles, the equation of inflaton is modifled 
as 

'cP + ^H^+im" + g^{x^))(t> = 0, (95) 
where the expectation value of x^ is defined as 

When the initial value of q is large as g > 3000, which corresponds io g> 3.0 x 10~^, the variance 
(x^) grows of order rr? j g^ and the backreaction onto the inflaton field cannot be ignored. This 
makes the oscillation of the inflaton incoherent, which flnally stops the resonance. The evolution 
of the variance (x^) for g = 3.0 x 10^^ is shown in fig. El 

The growth of the field perturbations can stimulate the enhancement of metric perturbations 
on sub-Hubble and even on super-Hubble scales |35i- Whether this enhancement occurs or not 
is sensitive to the evolution of the large-scale x fluctuation during inflation. In the model 
(|85|) . for example, strong amplification of the x fiuctuation requires a resonance parameter 
g ^ 1 at the beginning of preheating, in which case the large-scale Xk modes are suppressed 
during infiation !36i- Therefore it is difficult to amplify super-Hubble metric perturbations in 
this model. However there exist several models of inflation where metric perturbations can be 
amplifled during preheating j37| . This also provides an interesting possibility of overproducing 
primordial black holes in certain cases j38j . 
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Figure 4: The schematic diagram of the Mathieu chart and the typical paths for three types of 
resonance. The lined regions denote the instabihty bands (zeroth, first, second,- • •). Note that 
the field perturbation evolves in the region > 2q. 
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Figure 5: The evolution of the variance {x^) = {x^) /m? for g = 3.0 x 10~^. 

9 Summary 

In this lecture note 1 reviewed the concept of inflation, generation of density perturbations, and 
reheating after inflation. I explained how inflationary cosmology solves a number of cosmolog- 
ical problems such as flatness, horizon, and monopole problems. In addition inflation makes it 
possible to generate nearly scale-invariant density perturbations, which is consistent with obser- 
vations. Inflation is really an efficient mechanism to solve the cosmological problems associated 
with standard big-bang cosmology. In addition, elementary particles can be produced during 
the (p)reheating stage after inflation through the decay of the inflaton. 

So far there exist several cosmological scenarios that might be the alternative to inflation 
such as the pre-big-bang |39| and ekpyrotic/cyclic model IjiQj. Nevertheless, unlike the standard 
inflation, it is generally difficult to solve all of the major cosmological problems at the same time 
in these models without a fine tuning of the model parameters. It is fair to say that standard 
infiation with a slow-roll fiat potential is the most promising scenario in the very early universe 
among the models proposed so far. 

Nevertheless we still have unsolved problems even in the inflationary cosmology: What is 
the origin of the infiaton field? What is the state of the universe before infiation? The initial 
singularity can be avoided? 

Fortunately we begin to get the high-precision observational data, which is expected to 
reveal the nature of infiation in details in near future. From the theoretical side, there are a 
lot of extensive works to try to construct viable models of infiation based on superstring and 
supergravity theories. It is really encouraging that we are living in an exciting era where the 
early universe models can be strongly constrained from the upcoming observational data. 
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